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ABSTRACT 

The existence of an anti-correlation between the equivalent width (EW) of the narrow core of the iron Ka line and the luminosity 
of the continuum (i.e. the X-ray Baldwin effect) in type-I active galactic nuclei has been confirmed over the last years by several 
studies carried out with XMM-Newton, Chandra and Suzaku. However, so far no general consensus on the origin of this trend has 
been reached. Several works have proposed the decrease of the covering factor of the molecular torus with the luminosity (in the 
framework of the luminosity-dependent unification models) as a possible explanation for the X-ray Baldwin eff'ect. Using the fraction 
of obscured sources measured by recent X-ray and IR surveys as a proxy of the half-opening angle of the torus, and the recent Monte- 
Carlo simulations of the X-ray radiation reprocessed by a structure with a spherical-toroidal geometry by Ikeda et al. (2009) and 
Brightman & Nandra (201 1), we test the hypothesis that the X-ray Baldwin eff'ect is related to the decrease of the half-opening angle 
of the torus with the luminosity. Simulating the spectra of an unabsorbed population with a luminosity-dependent covering factor of 
the torus as predicted by recent X-ray surveys, we find that this mechanism is able to explain the observed X-ray Baldwin eff'ect. Fitting 
the simulated data with a log-linear L2-iokev - EW relation, we found that in the Seyfert regime (L2-iokev < 10"^^^ erg s~0 luminosity- 
dependent unification produces a slope consistent with the observations for average values of the equatorial column densities of the 
torus of XogN^ > 23.1, and can reproduce both the slope and the intercept for XogN^ ^ 23.2. Lower values of are obtained 
considering the decrease of the covering factor of the torus with the luminosity extrapolated from IR observations (22.9 < log < 
23). In the quasar regime (L2-iokev > 10^^^ erg s~^) a decrease of the covering factor of the torus with the luminosity slower than that 
observed in the Seyfert regime (as found by recent hard X-ray surveys) is able to reproduce the observations for 23 .2 < log < 24.2. 

Key words. Galaxies: Seyferts - X-rays: galaxies - Galaxies: active - Galaxies: nuclei 



1, Introduction (e.g. lYaqoob & Padmanabhanl2004l) . which points to the line be- 

, ^ . t . n iiig produced in cold neutral material. This material has often 

CO The observational signatures of reprocessed radiation in the X- identified as circumnuclear matter located at several thou- 

^ • ray spectra of Active Galactic Nuclei (AGN) are mainly two: a ^^^^ gravitational radii from the supermassive black h ole, and 

' Compton hump peaking around 30keV and a fluorescent iron ^^^1 ^^^^^^^ ^^^^ molecular torus (e.g., 

. ^ Ka line. While the Compton hump is produced only if the re- ,2006). although a c ontribution of the outer part of the disk 

K> processing matenal IS Compton-thick (CT, A^h > lO^^cm-^) (e.g. . iPetrucci et al.1 l2002h or of the broad-line region (BLR, 

q . the i ron Ka line is created also in Compton-thin matenal e.g., iBianchi et al."2008) cannot be excluded. In a recent study, 

_C^_. (e.g.,| Mattetal. || 2003|) . The iron Ka line is often observed as (20ll) found that the weighted mean of the ratio be- 

the superposition of two different components: a broad and a t^gg„ j^e FWHM of the narrow Fe Ka line and that of optical 

narrow one The broad component has a full width at half- ^^^^^ produced in the BLR is - 0.6. This implies that the size of 

maximum (FWHM) of > 30 000 km s'^, and is thought to Ka-emitting region is on average ~ 3 times that of the 

be created close t o the b ack hole in the accretion disk (e.g., gLR, and points towards most of the narrow iron Ka emission 

| Fabian etal. || 2000|) , or to be related to the presence of features produced in the putative molecular torus, 
created by partially covering warm absorbers m the line of sight 

(e.g., [Turner & Miller 2009, Mivakawa et aDl2012l) . While the One of the most interesting characteristics of the narrow 
broad component is observed in only ~ 35 - 4 5% of bright component of the iron Ka line is the inverse correlation be- 
nearby AGN (e.g., ' de La Calle Perez et al.l l2010h . the narrow twee n its equivalent width (^W ) and the X-ray luminosity 
(FWHM ~ 2. 000 km s"^ IShu et al.l |201o|) cor e of the iron (e.g., [Iwasawa & Taniguchj |1993|). The existence of an anti- 
line h as been foun d to be almost ubiquitous (e.g., [ Nandra et al.1 correlation between the equivalent width of a line and the 
120071 ISingh et al.l l201lh . This component peaks at 6.4 keV luminosity of the AGN con t inuum was found for the first 

time in the UV by iBaldwinI (1 19771) for the CIV .11549 line. 

Send offprint requests to: ricci@kusastro.kyoto-u.ac.jp and was then dubbed the Baldwin effect. A similar trend 
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Table 1. Summary of the most recent studies (along with the original work of llwasawa & Taniguchill 19931) of the X-ray Baldwin 
effect. The model commonly used to fit the data is logEW = A -\- 51ogLx,44, where Lx,44 is the 2-10 keV luminosity in units of 
10"^"^ erg s"^ The table is divided in works for which the fit was done per observation (i.e. using all the available observations for 
every source of the sample) and per source (i.e. averaging the values of EW and Lx of diflTerent observations of the same source). 
The table lists the (1) reference, (2) the value of the intercept and (3) of the slope, (4) the number and the type of objects in the 
sample, and (5) the observatory and instrument used. Dashes are reported when the value is not available. 



(1) 


(2) 


(3) 


(4) 


(5) 


Reference 


A 


B 


Sample 


Observatory/instrument 


Fits per Observation 


Shu et al. 2012 


1.58 + 0.03* 


-0.18 + 0.03 


32 (RQ) 


ChandrafHEG 


Shu et al. 2010 


1.58 + 0.03* 


-0.22 + 0.03 


33 (RQ+RL) 


Chandra/HEG 


Bianchi et al. 2007 


1.73 + 0.03 


-0.17 + 0.03 


157 (RQ) 


XMM-Newton/EFIC 


Jians et al. 2006 




-0.20 + 0.04 


101 (RL+RQ) 


XMM-Newton/EPlC + Chandra/UEG 


Jians et al. 2006 




-0.10 + 0.05 


75 (RQ) 


XMM-Newton/EFIC + Chandra/HEG 


Jimenez-Bailon et al. 2005 




-0.06 + 0.20 


38 (RQ) 


XMM-Newton/EFIC 


Zhou & Wans 2005 




-0.19 + 0.04 


66 (RQ+RL) 


XMM-Newton/EFIC 


Pase et al. 2004 




-0.17 + 0.08 


53 (RQ+RL) 


XMM-Newton/EFIC 


Iwasawa & Taniguchi 1993 




-0.20 + 0.03 


37 (RQ+RL) 


Ginga 


Fits per Source 


Shu et al. 2012 


1.64 + 0.03* 


-0.11+0.03 


32 (RQ) 


Chandra/HEG 


Shu et al. 2010 


1.63 + 0.04* 


-0.13 + 0.04 


33 (RQ+RL) 


Chandra/HEG 



Notes: * These values of the intercepts were corrected in order to be consistent with Eq.|5] as those reported in lShu et aT] ( 120101120121) 
are relative to the 2-lOkeV luminosity in units of 10^^ erg s~^ (^x,43)- 



was later found for several other emission lines such as Lyg , 
[C III] ^1908, Si IV ^139 6, Mg II ^2798 (.Dietrich et all l2002l). 
UV iron emission lines (iGreen et al mid-IR lines such 

as rARI IIU8.99A/m, [S IV] ^lO.Slyum and [Nell] ^12.8 lyum 
(iHonig et aL 2008), and forbidd en lines as [Oil] .13727 and 
[Ne V] ^3426 (ICroom et al.ll2002h . The slope of the Baldwin ef- 
fect for most of these lines has been shown to be st eeper for those 
originating from higher ionization species (e.g., iDietrich et al.l 
l2002h . The origin of the Baldwin effect is still unknown, al- 
though several possible explanations have been put forward, and 
it might be different for lines originating in different regions of 
the AGN. For the lines produced in the BLR it might be related 
to the lower ionization and photoelectric he ating in the BLR gas 
of more luminous objects (e.g. jNetzer et al.lll992b . 

In the X-rays, Iwasawa & Taniguchi (1993) using Ginga 
observations of 37 AGN, found the first evidence of an anti- 
correlation between the equivalent width of the iron Ka line and 
the 2 - lOkeV luminosity (EW oc L-^-20±o.04>^ j^^^ ^^^^^ 1^ 
ally called the X-ray Baldwin effect o r the Iwasawa -Taniguchi 
effect. Using A5'CA observations, N andra et al.l(ll997h confirmed 
the existence of such an anti-correlation, and argued that most 
of the effect could be explained by variations of the broad com- 
ponent of the iron Ka line with the luminosity. The advent of 
XMM-Newton, Chandra and Suzaku made however clear that 
most of the observed X-ray B aldwin effect is due to the narrow 
core of the iron Kg li ne (^e.g. JPage et"aDl2004 IShu et al.ll201QL 
iFukazawa et al .1120 111). The significanc e of the effect was ques- 
tioned by " Jimenez-B ailon et aP (|2005), who discussed the pos- 
sible importance of contamination from radio-loud (RL) AGN, 
which have on average larger X-ray lum inosities and smaller si 
natures of reprocessed radiation (e.g., iReeves & Tume 



However jGrandi et"an(l2QQ6b . studying B eppo SAX obso^ryations 
of radio-loud AGN, also found evidence of an X-ray Baldwin 
effect. This, together with the study of a large sample of radio- 
quiet (RQ) AGN carried out by Bian chi et all (7007) using 
XMM-Newton data , showed that the X-ray Baldwin effect is 
not a mere artifact. I Jiang et al] (l2006b suggested that the X-ray 
Baldwin effect might be due to the delay between the variabil- 
ity of the AGN primary energy source and that of the repro- 
cessing m aterial loca t ed far t her a way. Using Chandra/HEG ob- 
servations, IShu et al.l (l2010i 12012) confirmed the importance of 
variability, showing that averaging the values of EW obtained 
over multiple observations of individual sources would signif- 
icantly attenuate the anti-correlation from EW oc 



0.18±0.03 



to 



iilersig- 
^ 120001) . 



EW oc L-^-^^^^-^\ Chandra/HEG observations (IShu et al.ll2010h 
have also shown that the normalization of the X-ray Baldwin 
effect is lower (i.e. the average EW of the narrow Fe Ka compo- 
nent is smaller) with respect to that obtained by previous works 
performed using XMM-Newton. In Table[T]we report the values 
of the slope and the intercept obtained from the most recent 
XMM-Newton/EPlC and Chan dra/HEG works, along with the 
original Ginga work of Iwasawa & Tanig uchi (1993). 

An intriguing possibility is that the X-ray Baldwin effect 
is related to the decrease of the covering factor of the torus 
with luminosity (e.g.. Page et al. 2004, Zhou & Wang 2005), in 
the frame of the luminosity-dependent unification schemes. The 
equivalent width of the iron Ka line is in fact expect ed to be pro- 
portional to the coverin g factor of the torus (e.g., iKrolik etaP 
11994 llkeda et al.l 12009 ). so that a covering factor of the torus 
decreasing with the luminosity might in principle be able to 
explain the X-ray Baldwin effect. Luminosity-dependent unifi- 
cation models have been proposed to explain the decrease of 
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Fig. 1. Schematic representation of the torus geometry consid- 
ered. The angle 6i is the inclination of the line of sight with re- 
spect to the torus axis, while ^oa and are the half-opening 
angle and the equatorial column density of the molecular torus, 
respectively. H is the maximum height of the torus, while Rin 
and R are its inner and outer radius, respectively. 



Fig. 2. Variation of the fraction of obscured sources (/obs) and of 
the half-opening angle of the torus (^oa) with the 2-10 keV lu- 
minosity for^ of the most recent medium X-ray (2-l OkeV , 



iHasinged I2Q08I). hard X-rav (15-55 keV, Bu rlon etal.l l2Qllh 



and IR (iMaiolino et al.|[2007h studies (Eq.IUO 



the fraction of obscured objects (/obs) with the increase of the 
AGN output power. The first suggestion of the existence of a 
relat ion between fobs and the luminosity came about 30 years 
ago (iLawrence &Elvislll982l) . Since then the idea that the cov- 
ering factor of the obscuring material decreases with luminos- 
ity has been ga ining observationa l evidence from radio (e.g., 
'Grimes et al. 2004), infrared (e.g., Treister et al.'2008' "Mor et alj 
,2009, Gandhi etal. 2009), optical (e.g., Simpson 2005) and X- 
ravfe.g.. lUeda et al.l 120031 iBeckmann et alTl2009nUeda et al.l 
l201lh studies of AGN. Although luminosity-dependent unifica- 
tion models have been suspected for long to play a major role 
in the X-ray Baldwin eflfect, so far no quantitative estimation of 
this eflTect has been performed. 

In this work we quantify for the first time the influence of 
the decrease of the covering factor of the torus with the lumi- 
nosity on the equivalent width of the iron Ka line. Using the 
Monte-Carlo spectral simulations of a torus ^ yith a variab l e half- 
ope ning angle (^qa) recently pr esented by llkeda et all (l2009h 
and lBrightman & Nandral (1201 ll) . together with the most recent 
and comprehensive observations in diflTerent energy bands of the 
decrease of /obs with the luminosity, we show that this mecha- 
nism is able to explain the X-ray Baldwin eflTect. The paper is 
organized as follows. In Sect.[2] we present in detail our spec- 
tral simulations, and in Sect.[3]we compare the slopes and inter- 
cepts obtained by our simulations with those measured by high 
sensitivity Chandra/HEG observations, as they provide the best 
energy resolution available to date. In Sect. |4] we discuss our re- 
sults, and in Sect.[5]we present our conclusions. 



2. Simulations 

2.1. The relation between Oqa sind Lx 

As a proxy of the relationship between the half-opening an- 
gle of the molecular torus ^oa (see Fig.[T]) and the intrinsic X- 
ray luminosity of AGN we used the variation of the fraction 
of obscured sources /obs with luminosity measured by the re- 
cent medium (2- lO keV) and hard (15-5 5 keV) X-ray surveys of 
IHasinged (l2008h and lBurlon et"an (l2011l) . In the 2-10 keV band. 



from a combination of surveys performed by HEAO-1, ASCA, 
BeppoSAX, XMM-Newton and Chandra in the luminosity range 
42 < logLx < 46, Hasinger (2008) found 

/obs --0.2261ogLx + 10.342, (1) 

where Lx is the luminosity in the 2-10 keV energy range. In the 
15-55 keV band, using Swift/BAI to study AGN in the luminos- 
ity range 42 < logLx < 45, and fitting the data with 

/obs - ^low ^"^"^/^^ + /?Mgh (1 - ^-^"^/^^), (2) 

where Lrx is the luminosity in the 15-55 keV band. lBurlon et al.l 
(2011) obtained Riow = 0.8, R^gh = 0.2 and Lq = 
IQ43.7 g-i ^Ys Q ^gg^ ^i^g results of the IR work of 
iMaiolino et al.l (l2007l) . who found an anti-correlation between 
the ratio ^L^(6.7/im)/^L^(5 lOOA) and the [O III]/15007A line lu- 
minosity. This trend was interpreted as an eflTect of the decrease 
of the covering factor of the circumnuclear dust as a function of 
luminosity, with /obs varying as 

■^obs - 0.414' 
^ ^opt 

where Lopt = (ALa(5100A))/10'^^-^\ FoUowing IMaiolino et al.1 
(2007") it is possible to rewrite Lopt as a function of the X-ray 
luminosity: Lopt = L^^^^IO^^-^^. 

The fraction of obscured sources at a given luminosity can 
be easily related to the half-opening angle of the torus using 

^OA = COS-^/obs). (4) 

For Eq.|2] we converted the hard X-ray luminosity to the lumi- 
nosit y in the 2-10 keV band assuming a photon index F = 1.9 
(e.g., iBeckmann et ani2009l) . The three ^oa - Lx relationships 
used in this work are shown in Fig.O 

2.2. Spectral simulations and fitting 

llkeda et al.l (|^09) recently presented Monte-Carlo simulations 
of the reprocessed X-ray emission of an AGN surrounded by a 
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Energy (keV) 



Fig. 3. Top panel: extract in the 5.5-7.5 ke V region of a spectrum 
simulated using the model of Ikeda et al. I (l2009.) . The model has 
a torus with an half-opening angle of ^oa = 46. 2° (equivalent t o 
logLx = 42.7 according to the relationship of iHasingerl 12008 ), 
an equatorial column density of log A/^J = 23.8, and an inclina- 
tion angle of = 1°. The continuous line represents the fit to 
the simulated spectrum using for the continuum the same model 
we used for the simulations, and a Gaussian line for the iron K(2 
fluorescent line. Bottom panel: contribution to the chi- squared 
for the best fit to the simulations. 



^ — ^ — I ^ — ^ — I ' — ' — I 

100- 
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^ ^ I ^ ^ I , ^ I ^ ^ 

0.001 0.010 0.100 1.000 10.000 

2-10 keV X-ray luminosity (lO"*"* erg s"^) 

Fig. 4. Equivalent width of the iron Ka line versus the X-ray 
luminosity obtained simulating a torus with an equatorial col- 
umn density of logA^^ = 23.1 and a covering factor decreasing 
with the lum i nosity . The ^oa - Lx relationship us ed here is that 
of Hasinged (l2008 l) (see Fig. [2]). The relat i ons o f iBianchi et al.l 
( 2007) (black dashed line) and Shu et al.l (| 2012|; blue line, ob- 
tained averaging different observations of the same source) are 
also shown for comp arison (see also Tab le[T|). The normaliza- 
tion of the relation of iBianchi et al.l (l2007h has been fixed to an 
arbitrary value for comparison. 



three-dimensional spherical-toroidal structure. The simulations 
were performed using the ray-tracing method, taking into ac- 
count Compton down- scattering and absorption, and are stored 
in tables, so that they can be used to perform spectral fitting. The 
free parameters of this model are the half-opening angle of the 
torus ^oA, the line of sight inclination angle 6i, the torus equa- 
torial column density and the photon index T of the primary 
continuum. Note that should not be confused with the ob- 
served hydrogen column density Nu- In all objects, Nu < N^, 
with Nu = A/^H ^^^y ~ Fig.[T]a schematic repre- 

sentation of the geometry considered is shown. The dependence 
of Nu on the inclination angle is given by Eq. 3 of Ikeda et aO 
(12009). In the model of Ikeda e t al.l (12009 ) the ratio RiJR is fixed 
to 0.01, and the inclination angle of the observer 6i can vary be- 
tween 1° and 89°, A^^ between lO^^ and lO^^cm'^ while ^oa 
spans the range between 10° and 70°. Because of the assumed 
dependence of the half-opening angle on the luminosity, the in- 
terval of values of ^oa allowed by the model limits the range of 
luminosities we can probe. In the following we will consider lu- 
minosities above log Lx = 42, as below this value few AGN are 
detected, and several works poi nt towards a possible disappear- 
ance of the molecular torus (e.g. jElitzur & Sh losman 2006). The 
upper-limit luminosity we can reach depends on the ^oa - Lx re- 
lati onship used, and considering = 70 °, is log L max = 44.2 
for lHasinged (l2008h . lo g L^ax = 43.8 for iBurlon et al.l (l201lb 
and log L max = 45.0 for' Maiolino et al.l (2001). Due to its limi- 
tations in the values of ^oa permitted, the model of Ike da et al.l 
(l2009l) restricts most of the simulations to the Seyfert regime 
(logLx<logL^=44.2). 

To extend our study to the quasar regime (Lx > L^) we 
simulated the X-ray B aldwin eflTect using the spectral model of 
iBrightman & Nandra (201 1). This model considers the same ge- 
ometry as Ikeda et al. ( 2009), but has 7?in = 0, a -independent 
A^H, and it aflows diflTerent values of ^oa (26° to 84°). The 
higher maximum half-opening angle permitted by this model 



allows to reach higher maximum luminosities in the simu- 
lations for the relationships of Hasinger (2008; log L max = 
45. 3) and iMaiolino et al.l (|2007; logLmax = 46.2), while that 



of iBurlon et al.l (l201 11) is flat (/obs - 0.2) in the quasar 



regime. However, the lowe r boundary of ^oa in the model of 
Brightman & Nandra (2011) limits the lower luminosity we can 
reach in the simulations^ i n parti cular logLmin - 43 for the 
relat ion of [M aiolino e t aP (l2007h . Thus, we used the model 
of ilkeda et al.i (.2009.) to study the X-ray Baldwin eflTect in 
the Seyfert regime (42 ^ log Lx < 44.2), and the model 
of Brigh tman & Nandral (1201 ll) to probe the quasar regime 
(log Lx"> 44.2). 

To estimate the influence of the decreasing covering factor 
of the torus on the equivalent width of the iron line, we simu- 
lated, using the physical torus models, a large number of spectra 
using the ^oa - Lx relationships reported in Eq.[T]|3] Using the 
three diflTerent ^oa - Lx relationships, we extrapolated the value 
of ^oA for each luminosity bin (AlogLx = 0.01). We fixed the 
equatorial column density of the torus N^ to 26 diflTerent val- 
ues, spanning between 10^^-^ cm"^ and 10^^ cm"^, with a step 
of A log A/^ J = 0.1. To simulate an unabsorbed population, sim- 
ilarly to what is usually used to determine the X-ray Baldwin 
eflTect, for each value of Lx and N^ we considered inclination 
angles 6i between 1° and ^oa, with a binning of A^i = 3°. In 
XSPEC 12.7.1 (.Amaud 1996) we simulated spectra for each of 
the three ^oa - Lx relationships and for each bin of column den- 
sity. For the co ntinuum we used a power law with F = 1.9 (e.g., 
iBeckmann et al. 2009). Our choice of the photon index does 
not aflTect significantly the simulations, and adding a scatter of 
AF = 0.3 we found EW - Lx trends consistent with those ob- 
tained without scatter. The metallicity was set to the solar value 
in all our simulations, and the value of the normalization of the 
reflected component was fixed to that of the continuum. 

We fitted the simulated spectra in the 0.3-lOkeV band us- 
ing the same model adopted for the simulations, substituting the 
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22.5 23.0 23.5 24.0 24.5 25.0 22.5 23.0 23.5 24.0 24.5 25.0 

logNn logNn 

Fig. 5. Left panel: value of the slope {E) of the X-ray Baldwin effect obtained simulating the variation of the reprocessed X-ray 
radiation with the luminosity for tori with different values of the equatorial column density (in the range \ogN^ = 22.5 - 25), 
for the different ^oa - Lx relationships shown in Fig.[2](Eq.[T][3l) in the Seyfert regime (42 < log Lx < 44.2). The model used is that 
ofljkeda et al. ( 2009), and the simulated data were fitted with Eq.\5\ The re d lin e represents the value of the slope (B = -0.12 ± 0.04) 
obtained by fitting in the Seyfert regime the data per source reported in iShu et al.l(l2010k : the red shadowed area represents the Icr 
contour of the slope. Right panel: same as left panel, but considering the intercept (A) of the X-ray Baldwin effect obtained by our 
simulations. The red line and shadow represent the value of the intercept obtained by the fit to the Chandra/HEG data in the Seyfert 
regime and its Icr error (A = 1.64 ± 0.05), respectively. 



iron Ka line component with a Gaussian line. We obtained a 
good fit for all the simulations, with a reduced chi- squared of 
xl < 1.1. We report in Fig. [3] an example of a typical fit. The 
model used for the continuum fitting does not affect signifi- 
cantly the results, and iising a n alternative model like pexrav 
(iMagdziarz & Zdziarskilll995l) the equivalent width of the iron 
Ka line diflTers on average of only ~ 4%. We evaluated the 
equivalent width of the iron Ka line and studied its relation- 
ship with the luminosity. As an example, in Fig.|4lthe simulated 
X-ray Baldwin effect obtained using the model of llkeda et al.l 
(i2009) and the ^oa - Lx relationship of Hasingej (l2008h for 
logN^ = 23.1 is shown, together with the fit to the X-ray 
Baldwin ef fect o btained by the recent works of iBianchi et al.l 
(120071) andlShu et al. (2012). The spread in EW for a given lu- 
minosity is due to the range of values of Oi we considered: larger 
values of EW usually correspond to lower values of 6i. 

Most of the studies of the X-ray Baldwin effect performed in 
the last years have used a relationship of the type 

logEW = A-^B log Lx,44, (5) 

to fit the EW - Lx trend, where Lx,44 is the luminosity in units of 
10"^"^ erg s"^ . In order to compare the simulated the EW-Lx trend 
with that observed in unabsorbed populations of AGN, we fitted, 
for each value of N^, the simulated data with Eq.[5] using the 
weighted least-square method, with weights of w = sin^i. This 
allows to account for the non-uniform probability of randomly 
observing an AGN within a certain solid angle from the polar 
axis. 

3. The X-ray Baldwin effect 

Most of the studies on the X-ray Bald win effect have found, 
using Eq.[5l a slope of B - -0.2 (^e.g.. [iwasawa & Taniguchil 
I1993L lPageetal.1 [20041 iBianchi et al.l l2007h . However, in all 
these works the values of EW and Lx are obtained from in divid- 
ual ob servations of sourc es, and, as pointed ou t by Jiang et akl 
(l2006h and confirmed bv lShu et al.l (l2010l IMa) . flux variabil- 
ity might play an important role in the observed X-ray Baldwin 



eflTect. IShu et al.l (|2012|) recently found a clear anti-correlation 
between the equivalent width of the iron Ka line and the lumi- 
nosity for individual sources which had several Chandra/HEG 
observations. Using a sample of 32 radio-quiet AGN, they also 
found that the fit per source (i.e. averaging diflTerent observa- 
tions of the same object) results in a significantly flatter slope 
(B = -0.11 ± 0.03) than that done per observation (i.e. using 
all the available observations for every source of the sample; 
5 = -0.18 + 0.03). 

In order to compare the slope obtained by the simulations 
in the two diflTerent luminosity bands with real d ata, we fitted 
the data per source (33 AGN) of IShu et al.l (120101: obtained fix- 
ing cr = 1 eV) with Eq.[5] in the Seyfert and quasar regime. 
Consistently to what reported inl Shu et aP (12010*) we did not use 
the 3 AGN for which only upper limits of EW were obtained. 
Using the 28 AGN in the Seyfert regime we obtained that the 
best fit to the X-ray Baldwin eflTect is given by A = 1.64 ± 0.05 
and B = -0.12 ± 0.04. Chandra/HEG observations are available 
for only 5 objects in the quasar regime, and at these luminosities 
we obtained A = 1.5 ± 0.3 and B = -0.16 ± 0.22. 



3.1. Seyfert regime 

Fitting the simulated data with Eq.[5] in the Seyfert regime, we 
found that the slope obtained becomes flatter for increasing val- 
ues of the equatorial column density of the torus (Fig. 13 left 
panel). This is related to the fact that the iron Ka EW is tightly 
connected to both A^^ and ^oa, so that EW is more strongly de- 
pendent on ^oA (and thus on its variation) for low values of N^, 
which results in a steeper slope. On the other hand, when 
increases the dependence on ^qa becomes weaker, and the slope 
flatter. As shown in Fig. [5] (left panel), \ yhile the valiies of B ob- 
tained by the 6(^ a - Lx relationships of iBurlon et al.l (1201 ih and 
[Hasinger(2008) are similar along the whole range of consid- 
ered, flatter slopes are obtained for that of Maiolino et al. (20071). 
In particular for the latter relationship the correlation becomes 
positive for logN^ > 24. 
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Fig. 6. Left panel: value of the slope {B) of the X-ray Baldwin effect obtained simulating the variation of the reprocessed X- 
ray radiation with the luminosity for tori with different values of the equatorial column density A/^J, for the different ^qa - Lx 
relationships shown in Fig.[2](Eq.[T][3l) in the quasar regime (log Lx > 44.2). The model used is that of Bright man & Nandra[(l2Ql lb , 
and the simulated data were fitted with Eq.[5] The red line rep resents the value of the slope {B = -0.16 ± 0.22) obtained by fitting 
in the quasar regime the data per source reported in lShu et alT d^OlO): the red shadowed area represents the Icr contour of the slope. 
Right panel: same as left panel, but considering the intercept (A) of the X-ray Baldwin effect obtained by our simulations. The red 
line and shadow represent the value of the intercept obtained by the fit to the Chandra/HEG data in the quasar regime and its Icr 
error (A = 1.5 ± 0.3), respectively 



By combining the observations of the X-ray Baldwin effect 
with our simulations we can extrapolate the average value of the 
equatorial column density of th e torus of the unobscured AGN 
in the Chandra/HEG sample of IShu et"al ] ([2010). A slope con- 
sistent within Icr with our fit to the X-ray Baldwin effect in the 
Seyfert regime is obtained for > 10^^-^ cm"^ for the ^oa - Lx 
relationship of iHasingeJ (l2008l) . for A^^ > 10^^^ cm"^ for that of 
iBurlon et al (20 11), and for A^^ < 10^^ cm"^ for the relation- 



ship of Maio lino"et al.l (12007*). Comparing the value of the in- 
tercept obtained by the simulations with that resulting from our 
fit of the X-ray Baldwin effect in the Seyfert regime, we found 
that only a narrow range of average equatorial column densi- 
ties of the torus can reproduce the obser vations (right panel of 
Fig.O. For the ^oa - Lx relationship of iMaiolino et al.l (l2007h 
we found 22.9 < log N^ < 23, while for those of iBurlon et al.l 
(2011) and Hasinger (2008) we obtained logN^ ^ 23.2. Using 
both the values of A and B, it is possible to extract the average 
values of that can explain the X-ray Baldwin effect. As it 
can be seen from the tw o fig ures for bot h the ^oa - Lx relation- 
ship of Hasinger (l2008h and lBurlon~ al. ( 2011) the only value 
of column density consistent with both the observed intercept 
and slope is lo^N^ ^ wh ile for the IR ^oa - Lx rela- 

tionship of iMaiolino et al.l (l2007l) the values allowed are in the 
range 22.9 < logA^^ < 23. A similar result is obtained studying 
the A/5 chi- squared contour plot of our fit to the Chandra/HEG 
data. Using the model of Brightman & Nandral (l201ll) to simu- 
late the X-ray Baldwin effect in the Seyfert regime we obtained 
a range of A^^ consistent with that found using the model of 
llkeda et al.1 (5009) (23.2 < log A^^ < 233). The lower values of 
neede d to explain the X-ray Baldwin effect using the rela- 
tionship of IMaiolino et al.l (l2007l) are due to the larger values of 
/obs (and thus of ^oa) predicted by Eq.[3] This is again due to the 
fact that EW depends on both and ^oa, thus increasing the 
latter one would obtain lower values of the former. 



3.2. Quasar regime 

To study the behavior of B in the quasar regime, we fitted the 
simulated data using Eq.[5]for each value of for luminosities 

Lx > L^. Our simulations show that EW decreases more steeply 
in the quasar re gime than aUower l uminosities for the ^qa - Lx 
relationships of lHasingej(l2008h and lMaiolino et al.l(l2Q^ . with 
values of the slope ofB< -0.3 (left panel of Fig.[6l). The slopes 
obtained by the simulations are consistent with those found by 
fitting Chandra/HEG data in t he quasar regime for the ^oa - Lx 
relationship of Maiolino et al.l (12007) for A^ J > 10^^-^ cm"^. The 
values of the slope expected using the relationship of iHasingerl 
(2008) are steeper than the observed value for the whole range 
of column de nsities con s idered , while the flattening of the rela- 
tionship of IBurlon et al.l (l201ll) at high luminosities results in a 
slope of 5 ~ along the whole range of N^, consistent within 
Icr with the observations. The inter cepts obtained using the re- 
lationship of ^ Maiolino et al.l (l2007l) are consistent with the ob- 
served value for 22.5 < logN^ < 23.1 (right panel of Fig. [6]), 
with no overlap with the values of needed by the slope. We 
obtained intercepts that are consistent with the observations for 
22.7 < logA^J < 23.4 a nd 23.2 < logA^^ ^ ^4.2 for the relation 
of Hasinger (2008) and lBurlon et al.l (1201 1), respectively . Thus 
only the hard X-ray ^oa - Lx relationship Burlon et all (1201 lb 
is able to explain, for average values of the equatorial column 
density of the torus in the range 23.2 < logN^ < 24.2, at the 
same time both the intercept and the slope of the X-ray Baldwin 
effect at high luminosities. 

4. Discussion 

Since its discovery about 20 years ago, the existence of the X-ray 
Baldwin effect has been confirmed by several works performed 
with the highest spectral resolution available at X-r ays. So far, 
several possible explanations have been proposed. I Jiang et aP 
(2006) argued that the observed anti-correlation could be re- 
lated to the delay of the reprocessed radiation with respect to 
the primary continuum. The response of the circumnuclear ma- 
terial to the irradiated flux is not simultaneous, and one should 
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always take this effect into account when performing studies of 
reprocessed features as the iron Ka line or the Compton hump. 
[Sh u et al. ( 2010, 2012) have shown that averaging the values 
of Lx and EW for all the observations of each source results 
in a significantly flattened anti-correlation. However, by itself 
variability fails to fully account for the observed correlation 
(IShu et al.ll20T2h . 



4.1. Explaining the X-ray Baldwin effect with a 

luminosity-dependent covering factor of the torus 

A mechanism often invoked to explain the X-ray Baldwin ef- 
fect is the decrease of the covering factor of the molecular toru s 
with the luminosity (e.g.. Page et al. 2004, Bianchi et al."2007'). 
The decrease of the fraction of obscured sources with the lu- 
minosity has been reported by several works performed at dif- 
ferent wavelengths in the last decade, althoug h some discor- 
dant results have been presented (e.g., Dwelly ^ Pagel l2006l 
iLawrence & ElvislbOl OV In this work we have showed that the 
covering factor-luminosity relationships obtained in the medium 
and hard X-ray band can explain well the X-ray Baldwin eflTect in 
the 10"^^ - 10"^"^-^ erg s"^ luminosity range. In particular our simu- 
lations show that it is possible to reproduce the slope of the X-ray 
Baldwin eflfect with luminosity-dependent unification for aver- 
age equatorial column densities of the torus of logN^ > 23.1, 
and both the slope and the intercept for logN^ ^ 23.2 (Fig.O. 
In the same lii minosity range the ^oa - Lx IR relationship of 
iMaiolino et aP (12007 ) can explain the X-ray Baldwin eflfect for 
22.9 < log J < 23. 

In the quasar regime we ha ve shown that, wh ile the medium 
X-ray ^oa - Lx relationship of'HasingeJ (l2008h cannot explain 
the observations (Fig. [6]), the slope obtained by the IR ^oa - Lx 
relationship is consistent with our fit to the Chandra/HEG data 
in the same luminosity band. However, as for the latter relation- 
ship the range of values of required to explain the slope 
and the intercept do not overlap, it cannot be considered as a 
likely explanation. The hard X-ray ^oa - Lx relationship of 
iBurlon et^aP (l201ll) is flat above L^, and would thus produce 
a constant EW of the iron Ka line and a slope of 5 ~ for 
the whole range of N^, consistent within Icr with the value 
found using Chandra/HEG data. The intercept we found us- 
ing this relation is also consistent with the observational value 
for a large range of average equatorial column densities of the 
torus (23.2 < logA^^ < 24.2). The fraction of obscured AGN 
is not well constrained at high -luminosity, thus an y trend with 
the luminosity between that of Burlon et aP (12011*) and that of 
iHasinger [200 3) would be able to reproduce the observed slope. 
The relation of lHasingej(l2008h produces negative values of fobs 
for log Lx > 45.8, thus a flattening of the decline is expected be- 
low this luminosity. From Fig. 7 (left panel) in Hasinger (2008) 
one can see that indeed above log Lx - 45 the value of fobs ap- 
pears to be constant, similarly to what has been found at hard 
X-rays. A flattening of /obs in the quasar regime would also be 
expected when considering the large amount of accreting ma- 
terial needed to power the AGN at these luminosities. It must 
however be stressed that at high luminosities the X-ray Baldwin 
eflTect is not well studied, and the Chandra/HEG sample we used 
in this luminosity range is small, not allowing us to reach a firm 
conclusion on the variation of the iron Ka EW with the luminos- 
ity. Possible evidence of a flattening of the X-ray Baldwin eflTect 
in the quasar regime has been recently found by iKrumpe et al.l 
(I2010h . 



4.2. On the differences between X-ray and IR half-opening 
angle-luminosity relationships 

The main diflTerence between the diflTerent ^oa - Lx relationships 
we used is that those obtained in the X-rays are based on direct 
observati ons of the absorb ing material in the line of sight, while 
the one of lMaiolino et al.l (12007) is extrapolated from the ratio of 
the thermal infrared emission to the primary AGN continuum. 
The ^Q A - (5 100 A) relationship obtained by IMaiolino et al.l 
(l2QQ7h was converte d in ^qa - Lx using the L2kev - L/l(2500A) 
relation obtained by lSteflTen et al.l (l2006h . and then converted in 
the Lx - L^(5100A) relation assuming the optical-UV spectral 
slope obtained by Vanden Berk et al. ( 2001). All this is likely to 
introduce some error i n the ^qa - Lx obtain ed by their IR work. 
It has been argued by iMaiolino et al. l (2001) that the larger nor- 
malization of the ^oA - Lx relation they found is related to the 
fact that medium X-ray surveys such as that of Hasinged (l2008h 
are likely to miss a certain fraction of heavily obscured objects, 
which can instead be detected in the IR. This is due to the fact 
that for Compton-thick AGN at energies < lOkeV most of the 
X-ray emission i s depl eted. However, hard X-ray surveys as that 
of Burlon_etalJ (UpTT), which are much less biased by absorp- 
tion, have found a normalization of the ^oa - Lx relation consis- 
tent with that obtained at lower X-ray energies (see Fig.[T). Our 
results also show that the X-ray Baldwin efl'ect can be explained 
by a Compton-thin torus, and that larger values of /obs would 
imply even lower average values of N^, so that missing heav- 
ily obscured objects would not significantly aflfect our results. 
X-rays are also probably better suited to probe the material re- 
sponsible for the iron Ka line emission. The line can in fact be 
emitted by both gas and dust, and while the IR can probe only 
the latter. X-rays are able to infer the amount of both gas and 
dust. We thus conclude that for the purpose of our study the X- 
ray ^oa - Lx relations are better suited than those extrapolated 
from IR observations. 

It could be argued that the value of the hydrogen col- 
umn density commonly used to determine the fraction of ob- 
scured objects (and thus as an indicator of the presence of 
the torus) is A^h ~ 10^^ cm"^, while the torus is believed 
to have larger values of the equatorial column density, and 
that this obscuration might be related to the presence of dust 
lanes or molecul ar structures i n the h ost galaxy (e.g., iMattI 
pOOQ). However . iBianchi et al.l (200§), using the results of 
Delia Ceca et a D (I2008h . have showed that a similar decrease 
of fobs with the luminosity is obtained when setting this thresh- 
old to a larger value of line-of- sight column density. They found 
that the fraction of CT objects decreases with the luminosity as 
fcT ^ L"^-^^, similarly to what is found for fobs in the X-rays 
(e.g., Hasin ger 2008, see Eq.[T]). A similar result was obtained 
by lFiore et al.l (1200 9) from a Chandra and Spitzer study of AGN 
in the COSMOS field. This implies that the contribution of dust 
lanes or of galactic molecular structures to the observed Nw does 
not aflTect significantly the /obs - Lx relationships obtained in the 
X-rays. 

4.3. The equatorial column density of the torus 

The distribution of values of the equatorial column density of 
the torus in AGN is still poorly constrained. X-ray observa- 
tion can in fact infer solely the obscuration in the line of sight, 
and only studies of the reprocessed X-ray emission performed 
using physical torus model s lik e those of I keda et al . ( 2009|) , 
iBrightman & Nandral (1201 lb and lMurphv^ YaqoobI (12009) can 
help to deduce the value of the equatorial column density of the 
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Fig. 7. Iron Ka EW versus X-ray luminosities and predicted 
trends obtained for different values of the equatorial column den- 
sity of the torus. The points are the values of the equivalent width 
of the iron Ka line reported by Shu et al. ( 2010), obtained by av- 
eraging multiple Chandra/HEG observations of AGN. The two 
blue dash-dotted lines are the fits to our simulations of the X- 
ray Baldwin efifect using the ^oa - Lx relationship of lHasinge^ 
( 2008) for the Seyfert regime (log EW = 1.01 -0.17 logLx,44 for 
logA^^ ^ 22.5, and log £W = 2.05 -0.08 logLx,44 for log = 
25). The blue long-dashed lines represent the EW - Lx relations 
obtained in the quasar regime (logEW = 1.17 - 0.73 logLx,44 
for logA^^ ^ 22.5 and logEW = 2. 14 - 0.3 9 log L x,44 for 
logN^ = 25) using the relationship of Hasingei] (l2008l) . while 
the green d ashed lines r e presen t those obtained using the rela- 
tionship of iBurlon" et al.1 (I20T H). The intercepts obtained in the 
quasar regime have been modified in order to match those ob- 
tained at lower luminosities. 



torus. However, th i s kind of studies are still v ery scarce (e.g., 
iRivers et al.l I2OIII: iBrightman & Uedal l2012h . besides being 
largely geometry-dependent. Studies of AGN in t he mid-IR band 
perfor med using the clumpy torus formalism of iNenkova et al.l 
(l2008h have shown th at the number of clouds al ong the equator is 
A^o ~ 5 - 10 (iMor et al . 2009). If as reported bv lMor et al.l (|2009|) 
each cloud of the torus has an optical depth of ry ~ 30-100 
(i.e., logN^ ~ 22 - 23), the equatorial column density of the 
torus is expected to be logN^ ~ 22.5 - 24, in agreement with 
our results. It must however be noticed that the value of the in- 
tercept of the X-ray Baldwin eflTect obtained by the simulations 
is strongly dependent on our assumptions on the metallicity of 
the torus. This implies that our constraints of the average value 
of A^^ are also tightly related to the choice of the metallicity: 
lower values of the metallicity would lead to larger values of N^. 
To study this eflTect we repeated our study in the Seyfert regime 
using half- solar metallicities for the reflection model. We ob- 
tained that in this scenario, in order to explain the X-ray Baldwin 
eflTect, one needs values of the equatorial column density of 
the torus about two times larger than for the solar-metallic ity 
case (23.5 < logN^ < 23.7). To have a Compton-thick torus 
one would then need values of the metallicity of Z < 0.2 Zq. 
A Compton-thick torus has often been invoked to explain the 
Compton hump observed in th e spectrum of many unobscured 
AGN (^e.g. jBianchi et al.ll2004l) . It is still unclear however which 
fraction of the Compton hump is produced in the distant reflec- 
tor and which in the accretion flow. From our study we have 
found that in unobscured objects the X-ray Baldwin eflTect can 
be explained by a luminosity-dependent covering factor of the 



torus for an average value of the equatorial column density of 
logN^ ^ 23.2. This value is lower than the line-of-sight A^h of 
many Seyfert 2s (e.g., Ricci et al. 2011; Burlon et al. 2011), and 
it might be related either to the geometry we adopted or to pres- 
ence of objects with sub- solar metallicities. In particular, due to 
the constant reflection angle relative to a local normal in any 
point of the reflecting surface, the spherical-toroidal geometry 
produces larger values of EW (and thus lar ger values of the in- 
tercep t) with respect to a toroidal structure (iMurphy & YaqoobI 
I2OO9I) . 

It is possible that there exists a wide spread of equatorial 
column densities of the torus, thus one could envisage that this, 
together with the diflTerent values of Oi, would introduce the scat- 
ter observed in the anti-correlation. In Fig.|7]we show the X-ray 
Baldwin eflTect obtained from our simulations using the relation- 
ship of Hasinger (2008) for logN^ = 22.5 and logN^ = 25, in 
both the Seyfert and the quasar regime, to gether with the time- 
averaged Chandra/HEG data of Shu et all (l2010l) . From the fig- 
ure it is evident that all the data are well within the range ex- 
pected from our simulations, both in the Seyfert and in the quasar 
regime. It is important to remark that we do not know whether 
there exists a relation between the equatorial column density of 
the torus and the AGN luminosity. However, we have shown that 
the variation of the covering factor of the torus with the lumi- 
nosity alone can fully explain the observed trend, so that no ad- 
ditional luminosity-dependent physical parameter is needed. 

4.4. Luminosity-dependent unification of AGN 

The relation of fobs with the luminosity might be con- 
nected to the increase of the inner radius of the torus with 
the lumin osity due to du s t su blimation. Both near-IR rever- 
beration ( Suganumaet al. 2006) and mid-IR interferometric 
( Tristram & Schartmann 2011) studies have confirmed that the 
inner radius of the molecular torus increases with the luminosity 
as Rin oc L^-^. Considering the geometry of Fig.[T] the fraction 
of obscured objects (see Eq.|4]) would be related to ratio of the 
height to the inner radius of the torus (H/R) by 



^""^^ /? a|i + (////?)2- 

For values of H/R < 1, /obs oc H/R. In the orig inal formula- 
tion of the receding torus model (lLawrencdll99lb . the height H 
was considered to be constant. Assuming that R has the same 
luminosity-dependence as Rin, this would lead to /obs ^ L~^-^ . 
This has been shown to be inconsistent with the recent ob- 
servations, which point towards a flatter slope (/obs oc L"^-^^, 
e.g., Hasinger 20 08i se e Eq.[l] ), and which would imply H oc 
^0.25 Honig & BeckertI (1200 7*) have shown that in the frame of 
radiation- limited clumpy dust torus model, one would obtain 
H oc L^-^^, in agreement with the observations. 

The decrease of the covering factor with luminosity might 
also have important implications on the AGN dichotomy. It has 
been shown that Seyfert 2s appear to have on average lower 
luminosities and lower Eddington ratios (e.g., Beckmann et al] 
2009i lRicci et al.ll201 It) than Seyfert Is and Seyfert 1.5s, which 
suggests that they have on average a torus with a larger cover- 
ing factor. This idea is also supported by the fact that the lu- 
minosity function of type-II AGN has been fou nd to peak at 
lower luminosities than that of type-I AGN (Delia Ceca et al.l 
2008i iBurlon et al.l l201l|). and by the results obtaine d by the 
recent mid-IR work of Ikamos Almeida et al.l (l201ll see also 
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lElitzui]l2012h . A difference in the torus covering factor distri- 
bution between different types of AGN would also explain why 
the average hard X-ray spectrum of Compton-thin Seyfert2s 
shows a larg er reflection comp onent than that of Seyfert Is and 
Seyfert 1.5s (iRicci et al.ll201l]) . A similar result was found by 
iBrightman & Uedal (l2()12h : studying high-redshift AGN in the 
Chandra Deep Field South they found that more obscured ob- 
jects appear to have tori with larger covering factors, although 
they did not find a clear luminosity dependence. 



5. Summary and conclusions 

In this work we have studied the hypothesis that the X-ray 
Baldwin eflTect is related to the decrease of the covering factor 
of the torus with the luminosity. We have used the physical torus 
models oflk eda et al. (2009) and lBrightmaii& Nandra (201 1 ) to 
account for the reprocessed X-ray radiation, and the values of the 
fraction of obscured sources obtained by recent surveys in the X- 
rays and in the IR as a proxy of the covering factor of the torus. 
Our simulations show that the variation of the covering factor 
of the torus with the luminosity can explain the X-ray Baldwin 
eflPect. In the Seyfert regime (Lx < 10^^-^ erg s"^), the observed 
EW - Lx trend can be exactly (both in slope and intercept) re- 
produced by an average value of the equatorial column density 
of the torus of logN^ ^ 23.2, while a slope consistent with the 
observations is obtained for a larger range of column densities 
(logA^^ > 23.1). At higher luminosities (Lx > 10"^"^-^ erg s"^) the 
situation is less clear due to the small number of high-quality 
observations available. Moreover, it is not clear whether in the 
quasar regime /obs still decreases with the luminosity similarly 
to wh at is found in the Seyfert r egime (as shown b y 'Hasinger* 
l2008h . or is constant (as found bv lBurlon et~ani201lh . A flatten- 
ing of the ^oA - Lx relationship would be able to explain the 
observations at high luminosities (for 23.2 < logN^ < 24.2), 
and it might naturally arise from the large amount of accreting 

mass needed to power these luminous quasa rs. 

In the next years ASTRO-H (Takahashi e tani2010h . with its 
high energy-resolution calorimeter SXS, will allow to constrain 
the origin of the narrow component of the iron Ka line, being 
able to separate even better than Chandra/HEG the narrow core 
coming from the torus from the flux emitted closer to the central 
engine. ASTRO-H will also be able to probe the narrow iron Ka 
line in the quasar regime, allowing to understand the behavior of 
the X-ray Baldwin effect at high luminosities. 
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